Landau-Zener Tunnelling in Waveguide Arrays 
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Landau-Zener tunnelling is discussed in connection with optical waveguide arrays. Light injected 
in a specific band of the Bloch spectrum in the propagation constant can be transmitted to another 
band, changing its physical properties. This is achieved using two waveguide arrays with different 
refractive indices, which amounts to consider a Schrodinger equation in a periodic potential with a 
step. The step causes wave "acceleration" and thus induces Landau-Zener tunnelling. The region 
of physical parameters where this phenomenon can occur is analytically determined and a realistic 
experimental setup is suggested. Its application could allow the realization of light filters. 

PACS numbers: 42.82.Et; 42.25.-p; 03.65.-w 

When a quantum system is subject to an external 
force, a non-adiabatic crossing of energy levels can oc- 
cur. This phenomenon is known as Landau-Zener tun- Upper Band 
nelling P, Q and some of its recent observations are for 
Josephson junctions and optical effective two-level 
systems Q. On the other hand, the problem of quan- 
tum motion in a periodic potential was solved already 
in the 1920's (see e.g. and gives rise to band spectra 
and Bloch states. Nowadays, the observation of Landau- 
Zener tunnelling between Bloch waves is at the fronteers 
of research in Bose-Einstein condensates (BEG) in op- 
tical lattices H H H E The external forcing, 
responsible for Landau-Zener tunnelling, is created by 
either placing the BEG in a gravitational potential |^ or 
accelerating the optical lattice itself . 

In this Letter we propose a new way of generating 
Landau-Zener tunnelling. We consider waveguide ar- 
rays J where the periodic potential of the Schrodinger 
equation is provided by the spatial oscillation of the re- 
fractive index in the transversal direction. Tunnelling is 
caused by combining two waveguide arrays with differ- 
ent refractive indices (see Fig^. As we will see, this 
corresponds to creating a step in a periodic potential. 

For arrays of coupled waveguides, the longitudinal di- 
rection z (see Fig^), along which the refractive index 
is constant, plays the role of "time" in the station- 
ary regime. The refractive index varies only along the 
transversal direction x, which represents space in a 1 -I- 1 
(space- "time" ) dimensional picture. Various linear and 
nonlinear phenomena have been observed in waveguide 
arrays: discrete spatial optical solitons diffraction 
management 0| , excitation of Bloch modes , gener- 
ation of multiband optical breathers jl^ and of single 
band-gap solitons [l3> anomalous band-gap transmis- 
sion regimes [l^ . Fast progress in discovering various 
nonlinear effects in waveguide arrays has been possible 
due to the introduction of the tight-binding approxima- 



Lower Band 



FIG. 1: Schematic picture of the combination of two waveg- 
uide arrays of the same spatial period but with different re- 
fractive indices (shown by the different grey levels). Choosing 
appropriately [however, see formula below] the difference 
in refractive indices the intensity in the lower band mode is 
almost completely transferred to the upper band mode due 
to Landau-Zener tunnelling. 



tion [19II20II2H . which reduces the nonlinear Schrodinger 
equation to the discrete nonlinear Schrodinger equa- 
tion [13. 

However, such a reduction eliminates the rich band 
structure of the periodic medium and only a single Bloch 
band is left. On the contrary, we want to maintain the 
band structure and, hence, study transitions between the 
bands. Indeed, as we will see below, the coupling of two 
waveguides with different refractive indices introduces a 
step in the periodic potential and, consequently, Landau- 
Zener tunnelling. In the following we will study only 
transitions between the first two bands, denoting them 
upper and lower band. In particular, we propose to inject 
light into the left waveguide array with a given angle in 
order to populate the lower Bloch band, and retrieve it 
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FIG. 2; Approximation of the periodic potential with a step 
used in the Schrodinger equation describing the waveguide 
arrays with different refractive indices. 



from the right waveguide array. Below, we derive analyt- 
ically the step size bounds inside which most of the inten- 
sity of the lower band mode is transferred to the upper 
band mode, creating also a spatial separation between 
lower and upper band light. We demonstrate this effect 
by performing numerical simulations of the Schrodinger 
equation in a periodic potential with a step. 

The waveguide array refractive index profile in the 
transversal direction x, which has a periodic rectangular 
shape, is approximated by a harmonic potential, while 
the step in the refractive index is substituted by a slope 
—a of height A (see Fig. The functional form of the 
harmonic potential is selected in order to keep theground 
state at zero energy, as done in the case of BEC IT] . Such 
an approximation allows a simple analytical treatment of 
the problem. Thus, in the linear regime, the adimension- 
alized Schrodinger equation of the optical system can be 
written as follows (see e.g. Ref. 22] ) 

+ +(n.) + 2.sin^.)*^0, (1) 

where Vl/ stands for the complex envelope of the electric 
field and w is the height of the harmonic potential. More- 
over, 

V{x) = for a; < 0; V{x) = -A for x > A/a 
and V{x) = -ax for < a; < A/a. (2) 

Via the simple transformation "if 'i/ exp[izV{x)], 
wave equation gets a well known form (see Refs. 

mniiii): 
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where a plays the role of acceleration in the Landau- 
Zener phenomenon. However, we should keep in mind 
that acceleration takes place only within the step region 
< X < A/ a, unlike the previously considered cases of 
EEC's, where the whole condensate is accelerated (see 
e.g. Ref. Using a two mode approximation the 

wave function ^ is written as follows 



^ = a{z)e 



Kx 



h(z)e 



i(K-2)x 



(4) 



With our conventions the zone-boundary mode 
wavenumber is 1. It should be mentioned that the 
two mode approximation works better just in the 
vicinity of zone boundaries {K —^ 1), exactly where 
Landau-Zener tunnelling takes place. Following Ref. j^, 
we substitute Expr. Q into the wave Eq. |(2Jl. Assuming 
K — 1 and removing a common phase dependence in 
a{z) and b{z), we get the Landau-Zener model 0, Q in 
its original form 
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(5) 



Thus, according to Landau-Zener's result, tunnelling 
from the lower zone-boundary mode to the upper band 
takes place with the following rate 
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(6) 



We consider cases in which the acceleration constant a 
is large. Thus, the tunnelling rate is close to one, mean- 
ing that almost all light intensity approaching the step 
is transferred to the upper band. This is confirmed by 
numerical simulations. 

Outside the step, where acceleration is absent, one can 
write down the wave- functions and the dispersion rela- 
tions in simple approximate form (see e.g. [Sj). 
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(7) 
(8) 



where k = K — 1 is the wavenumber detuning from the 
zone-boundary and the -|-(— ) sign indicates the upper 
(lower) band. /3 is the dimensionless propagation con- 
stant. The dispersion relations © for the two bands 
are schematically shown in Fig.|21 The picture is similar 
to the one observed in experiments |l7j . Let us remark 
that at the zone-boundary (k — > 0), the amplitudes of 
the upper and lower band modes are — 2 sin a; and 
\'^-\ = 2cosa;, respectively. This means that the light 
intensity in the lower band is concentrated inbetween the 
waveguides centers, while, in the upper band, intensity is 
concentrated on waveguides centers. This property is a 
clear experimental indication whether the wave is in the 
lower or upper band. 

The experimental setup could be as follows. One 
should inject a lower band mode with non zero but small 
relative wavenumber k. This is accomplished by choos- 
ing for the light beam a direction forming an angle 6 with 
respect to the z direction such that tan 9 = k. Hence, the 
wave front will move towards the step. An analysis of the 
dependence of the tunnelling mechanism on the physical 
parameters appearing in Fig. 13 shows that the transition 



3 



t 




K[7i/period] " ^ ^ 

1 2 



FIG. 3; Schematic band-gap structure and picture of the 
Landau-Zener tunnelhng process, w is the gap between the 
bands (the height of the period potential), d is the width of 
the upper band and k is the initial detuning of the lower 
band mode wavenumber from zone-boundary. P is the di- 
mensionless propagation constant (see formula |SJ| . A is the 
height of the step. Initially, light is injected in the left ar- 
ray, populating the lower band mode(a). Going across the 
step K decreases, reaching zero as the mode approaches the 
zone-boundary causing Landau-Zener tunnelling (b). After 
tunnelling, most of the light intensity is transferred to the 
upper band mode, whose wavenumber decreases until the end 
of the step is reached. This is the light we observe in the right 
array (c). A smaller light intensity remains in the lower band 
and is observed in the left array (d) . 

to the upper band mode verifies only if the refraction 
index step A fulfills the following inequalities 

A/3 + w < A < Al3 + w + d, (9) 

where A/3 = /3-(0) — (3-{k) is the variation of the 
propagation constant between the initial state and zone- 
boundary and d is the width of the upper band. 

Let us try to justify this result by commenting at the 
same time the results of some numerical simulations. 
These are performed by fixing w — 0.5, k — 0.2 and vary- 
ing the step size A. Waveguide centers are placed every 
period, with the first waveguide at half a period from 
the left boundary. The refractive index step is placed in 
the middle of the array. If the refractive index step A is 
within the above limits @, the lower band mode cannot 
overcome the step and when it reaches the zone-boundary 
Landau-Zener tunnelling to the upper band occurs. Light 
is partially transmitted to the right in the upper band 
and reflected to the left in the lower band. This regime 
is demonstrated in Fig.^ which evidentiates, even visu- 
ally, the fact that light in the left array is concentrated 
inbetween waveguides, while it lies at waveguide centers 
on the right array [see also the form of wavefunctions 
Q]. If the step size is higher than the upper bound of 
Exp. the lower band mode cannot overcome the step 




FIG. 4: Landau-Zener tunnelling from the lower to the upper 
band. Step size A = 0.7 is taken within the limits of Exp. 
@. Step is placed at a; = 5 and lower band mode is injected 
into the first five periods. 

and total reflection takes place (see the upper graph of 
Fig. [SJ. On the other hand, if 

A<A(3 (10) 

the lower band mode is able to overcome the step and 
to penetrate to the right side without tunnelling (see the 
lower graph of Fig. |3J). Indeed, wave intensity is now 
concentrated in-between the waveguides and one can con- 
clude that only lower band modes are present in the array. 
Finally, if the step is located within the following limits 

Af] < A< Ap + w, (11) 

the wave on the left array has no counterpart on the left 
with the same propagation constant, thus, no stationary 
penetration of the light through the step is possible, more 
or less like in the upper Fig. \^ 

Concluding, a novel type of linear optical tunnelling 
effect is discovered. This is discussed in connection 
with waveguide arrays, which have a spatially oscillat- 
ing refractive index. We explain this effect resorting 
to Landau-Zener model, which is commonly used for 
accelerated quantum-mechanical two-level systems. In 
our case, tunnelling between different bands takes place 
while a wave passes through a step in the refractive in- 
dex. In waveguide arrays, Bloch bands in the propa- 
gation constant (light wavenumber) play the role of en- 
ergy bands in quantum mechanics. Numerical simula- 
tions show that, if certain limits in the refractive index 
step are respected, a spatial separation of light in differ- 
ent bands can be achieved. More interesting for appli- 
cations is, perhaps, the use of this mechanism to build 
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FIG. 5: Upper graph: Total reflection from the step. Step 
size A — 1 is taken higher than the upper bound of Exp. 
Q. Lower graph: Penetration through the step without tun- 
nelling. Step size A = 0.07 is taken within the step values 
form the inequality 1101 . 

light wavenumber filters. Indeed, by appropriately choos- 
ing the physical parameters, it could be possible to shift 
the central wavenumber of a light beam and to reduce 
the spread in both frequency or wavenumber. The inclu- 
sion of a weak nonlinearity does not qualitatively alter 
the picture, while strongly nonlinear cases need a sepa- 
rate treatment. In particular, the scattering process of 
optical gap solitons with the step could be of special in- 
terest. Moreover, it should be observed that the effect 
discussed in this Letter is generic for systems with pe- 
riodic potential and that applications in different fields 
could be found. For instance, the analysis of a similar 
process in Bose-Einstein condensates is in progress. 
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